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In Class. Quantum Grav. 26 (2009) 175006 Bizori et al discuss the power-law tail in 
the long-time evolution of a spherically symmetric self-gravitating massless scalar field in 
odd spatial dimensions. They derive explicit expressions for the leading order asymptotics 
for solutions with small initial data by using formal series expansions. Unfortunately, this 
] approach misses an interesting observation that the actual decay rate is a product of asymp- 

■ totic cancellations occurring due to a special structure of the nonlinear terms. Here, we 

kjpj! show that one can calculate the leading asymptotics more directly by recognizing the special 

^ ■ structure and cancellations already on the level of the wave equation. 
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Since the works of John []J and Asakura |2] who studied nonlinear wave equations with power 
nonlinearities it is known that the late-time asymptotics can be ruled by the nonlinearity even for 
solutions starting from small initial data (for which one might naively want to ignore the nonlinear 
terms as causing only "higher order" corrections). For wave equations with nonlinearities containing 
first derivatives the asymptotics may even depend on the particular linear combination of the terms. 
(3jT)! Christodoulou [3J and Klainerman [4| discovored special null structures, which lead to a faster than 

generic decay or decide about a global existence of solutions. The underlying mechanism is based 
on asymptotic cancellations of the leading order terms in these special nonlinear structures (cf. 0] 
for a detailed analysis of such cancellations). The purpose of this Comment is to demonstrate that 



the same phenomenon occurs here, in the wave equation for the scalar field, and make use of it to 
simplify the calculation of the leading asymptotics at late-times. 

In the commented Article 0], later referred to as BCR, evolution of a self-gravitating real 
massless scalar field is considered. The Einstein equations for a d + 1-dimensional metric with 
odd d > 3 restricted to spherical symmetry 
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ds 2 = e 2a(t,r) (jP&r)^ + ^ + (1) 



are analyzed, where dQ 2 i _ 1 is the round metric on the unit (d — 1) -dimensional sphere. The scalar 
field satisfies a (quasilinear) wave equation with smooth, and compactly supported initial data 
(0, 0)t=o = (sf,eg) where e is a small number. The functions 0,/3 and m = (1 — e~ 2a )r d ~ 2 are 
formally expanded in the Taylor series in e about Minkowski spacetime thq = /?o = </>o = and 
substituted into the field equations. This gives an infinite hierarchy of equations on 0j,/3j,mj, 
i = 1, 2, 3, ... which can be solved recursively. 
At the first order mi = (5\ = while 

□0i = O, (01, 0i)t=o = (/,<?) (2) 

can be solved explicitly. At the second order 02 = and 

m' 2 = Kr d - 1 {4>l + 4>' 2 ?) , (3) 
m 2 = 2Kr rf - 1 0i0' 1 , (4) 

A = (5) 
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where n = — . And at the third order 

d-l 

Ofa = 2fafa + fafa + (3' 2 fa. (6) 

In BCR it is claimed that the late-time asymptotics of <f> is dominated by that of fa. In order 
to calculate it some of the functions fa, fa, mi must additionally be expanded in powers of r _1 
and then inserted into ([5]). Here, cancellations of the leading terms occur since this equation has 
a very special structure of the nonlinear terms originating from the wave equation for 0. Hence, 
next-to-leading order terms become important in the asymptotics. 

Below, we present a method which makes use of these cancellations already at the level of 
the wave equation (cancellation of terms before their evaluation) and thus reduces the amount of 
necessary asymptotic information about the source functions on the right-hand side. In Section[I]we 
regroup the nonlinear terms, eliminate the subdominant ones and calculate the leading asymptotics 
for </>3 by solely evaluating the dominant nonlinear term. In Section [H] we prove estimates which 
give a rigorous background for the term selection in Section [H 



I. ANALYSIS OF THE LEADING ORDER ASYMPTOTICS 

A. 3+1 dimensions 

For d = 3 the solution of the free wave equation ([2]) can be written as (cf. (13)-(14) of BCR 
with I = 0) 

a(u) - g(y) 

fa(t,r) = . (7) 

r 

where u = t — r,v = t + r and the function a is determined by / and g and has compact support 
in [-R, +R\. The most obvious way of calculating the asymptotics of is to substitute the above 
function into (IS])-©, calculate fa, fa, fa, fa[, fa, 02 an d insert into ([6]) to obtain the desired decay 
in time for fa, as was done in Section III of BCR. However, as follows from the rough estimates 
(cf. those obtained in Section [H]) , the right-hand side of © is a function supported in the vicinity 
of the lightcone t = r and decaying like 1/r 3 . It suggests that fa should decay in time like 1/t 2 . 
As the section III of BCR shows, the true decay is by one power faster. Indeed, there happens a 
cancellation of leading terms in the asymptotic expansion. Here, we want to explain the cancellation 
mechanism already on the level of the differential equation by regrouping terms to form special 
structures. This transformation also allows for a considerable simplification of the calculations. 

By introducing null derivatives d± := \{dt ±<9 r ) and rearranging terms the equation © can be 
rewritten as 

Ufa = --fad-fa + 2d + (3 2 d + fa + -8-(r fad-fa) + -fad%{rfa). (8) 

The last term is identically zero for t > R because there fa is a purely outgoing wave ([7]). The 
second last term, in the process of inversion of □ = ^5_9+(r-), will turn out to be a complete 
derivative in the ingoing direction u. Since r(3 2 d-fa has compact support in u this term will 
vanish after integration. The second term of the above expression, as is explained in Section HU 
has faster decay in r (by at least one power) than any other combination of null derivatives and 
together with the compact support in u (localization near to the lightcone) leads to a faster decay 
in time for fa. Hence, it is the first term that will determine the asymptotic behavior of fa at late 
times (see Section |H] for proof). We denote it symbolically 

a fa * --fad- fa. (9) 
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Hence, all we need to calculate is the leading order behavior of /?2- Substituting ([7]) into ([3]) and 
integrating from r = we first find for t > R 



rri2(t,r) = 4tt 2 / a (x)dx 



a 2 (t - r) 



%irF{t-r) + 0{r- 1 ), 



(10) 



t-r 



where F{u) := Ii(u) is defined in (28) of BCR. For t > r + R we have exactly m2(t,r) = 0. Next, 
integrating equation from r = and using -F(it) = for u > R we get for r + R > t > R 



&M=87r / 7— [F^O + OC*" 1 )]^ 



"t 2 



(t - u') 2 

y F(n') d«' + 0(r 3 ) = ^G{t -r) + 0{t- ?J ) 



t-r 
oo 



(11) 



t-r 



where G(u) := ll(u). Otherwise, for t > r + R, we have P2{t,f) = 0. Substituting this and (J7]) 
into © and using (21) of BCR (with I = 0) we get 



+R t+r 



fo(t,r) = J ' drj J 



-R t-r 



G( V )a'(ri)+0 



(£ + ri) 



(12) 



Now, for t±r^>R elementary integration over £ and by parts over r/ yields the asymptotic result 
(34) of BCR 



fo(t,r) 



( t 2 _ r 2)2 



r + o 



t 2 — r 2 



T := -2 5 vr / F{u)a{u)du. (13) 



B. Higher dimensions 



At present we cannot rigorously prove the leading order asymptotics due to a lack of an optimal 
decay estimate for the wave equation in d+1 dimensions, but analogously to the the 3+1 case, 
we are able to regroup and estimate the right-hand side terms of the wave equation for ^3 and so 
determine the leading order source term (see Section [TTJ) . 

We can rearrange the nonlinear terms in the equation © to write them as 

□03 = ^-,3- (r^fca-fc) + (W&a, 0x) (14) 

The second source term of the above expression, as is explained in Section [TU has faster decay in 
r (by at least one power) than any other combination of null derivatives and together with the 
compact support in u leads to a faster decay in time for $3. Hence, it is the first term that will 
determine the asymptotic behavior of ^3 at late times. See Section |H] for a quantitative analysis. 
We denote it symbolically 



□03 = ^-,d- (r^&fl-fc) 



(15) 
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Using (21) of BCR we get 



1 



t+r r+oo 



t—r J — oo 



(e-??) 2/+2 /W 



(16) 



where I = {d — 3)/2. The inner integral over r/ can be integrated by parts (it produces no boundary 
terms since the integrand has compact support in rj) 



1 



2i+2 r ;+i 



t+r r+oo 

d£ / dr] d 



t—r J — oo 



i+i 



{t-ri? W fod v 4>X- (17) 



Now we only need to find the asymptotic form of the expression fad-fa. According to (13)-(14) 
of BCR we have 



fa(t,r) 



and 



' +0(r- 2 - 1 ), fa(t,r)<* t^t 1 +0(r- 2 - 1 ), fa(t,r)^- 



„i+/ 



-l+l 



+ 0(r 



-2-l\ 



oo 

m 2 {t,r)% 2k J [a^Xx)] 2 dx + Otr- 1 ) =: 2/cFj(t - r) + 0(r 

t-r 



(18) 



(19) 



where Fi(u) := ih^iu). Again, for t > r + R we have exactly m 2 (t,r) = 0. Next, integrating 
equation ([5]) from r = and using Fi (u) = for u > R we get for r + R > t > R 



fa(t,r) ^(d-2)2 K J (t _^ )2+2 , [F,(n / ) + 0(t- 1 )]^ 

t-r 

oo 



, (21 + 1) 16tt 
''(21 + 



(20) 



where Gi(u) := ll +1 (u), otherwise, for i > r + i?, we have fa(t,r) = 0. Substituting this and ([18 
into (|17p we obtain 



fa(t,r) 



(21 + 1) 2 2l+5 ir f t+r „ f + °° , n 
d£ dr] d n 



(21 + 2) r l + l 
Analogously to (22) or BCR it holds 

t+r 

f d£d n 



t-r 



l+l 



it + ^l ^^Hv). (21) 



Pl(») 



t-r 



l+l 



(Z + v) 



2+21 



,^ (21 + 2) r l +i 
{ ' (21 + 1) £3+3* 



1 + 



(22) 



which applied to the above integral yields the asymptotic result (44) of BCR 



h(t,r) = -^- + O f ' 



+oo 



r ; := (-iy +1 2 3 '+ 5 vr / FMa^^dr], (23) 



£3Z+3 ^ \^3Z+4 

for late times (t > r), where in Ti we integrated by parts over r/. The asymptotic formula (45) of 
BCR can be straightforwardly derived from (I2ip as well, it only requires some other asymptotic 
expansion of the integral (f22~|) . 
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II. ESTIMATES 
A. 3+1 dimensions 

In this section we will prove decay estimates for 1712, @2 and <f>3- They will lose information about 
the exact amplitudes but will be helpful in separating the leading asymptotics from the subleading 
corrections decaying faster. 

From d?l) we know that is supported in the strip < t — R < r < t + R and can estimate 
there its derivatives 

\Mt,r)\Mt,r)\,\<l>[(t,r)\<± |0 1 + ^|< J_ (24) 

where (x) := 1 + \x\ and "<" means "less or equal than" up to some multiplicative constant which 
we skip for brevity. Observe that for functions supported in the strip \t—r\ < R estimates by powers 
of (r) and (t) are equivalent, since there exist constants Ci,C 2 such that C\(t) < (r) < C 2 {t). 
From ([3|) we have m 2 (t, r) = for r < t — R and for < t — R<r<t + Rwe estimate 



\m 2 (t,r)\ < / r' z M^j 2 + |^| 2 *'< / —-dr' < C (25) 
JO v 7 Jt-R V I 

where C depends only on R. Actually, there is a universal bound \rri2(t,r)\ < M where M is the 
total energy ("ADM mass") of <j>\. Moreover, for small r, say r < 1, we immediately see that 
\ m 2{t,r)\ < r 3 . These two facts can be put together to give 

\m 2 (t,r)\<^ (26) 



Next, from ([5]) we get /3 2 (t, r) = for r <t — R and for < t — R < r < t + R we again estimate 
.„ , f r m 2 (t,r') f t+R 1 ,, 2R . 1 

We will also need a similar estimate for the outgoing derivative of [5 2 - Therefore we combine the 
equations (IS)-© to get the identity 

d+p 2 = = J o (</>i + 01) ^' - 2 J ^dr'. (28) 

Using the above estimates on (pi + §\ and m2 we find 

^ < ~Lw dr ' + Lw dr '^w (29) 

for < t - R < r < t + R as well as d + (3 2 (t, r) = for r < t - R. 

Finally, we analyze the various source terms in the wave equation ([U]). Let us split the solution 
into four components introduced in (JSj) 

D<f> 3A = --fod-fa, Ufo B = (30) 

r 

□0 3C = ^-(r/3 2 9_0 1 ), Q/> 3 d = ^/?2^N>i), (31) 
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with 03 = 03,4 + 4>3B + 03C + 03.D- All terms on the right hand side are supported in \t — r| < R. 
By the above bounds we can estimate the first two components 



-fad-<fti 



< 



1 



< 



1 



r(r)(t) 2 ~ r(r) 3 ' 



\d + p 2 d+(j)i 



< 



1 



< 



1 



(t) 3 (r) 2 ~ (r) 5 



Now, e.g. from 0], we find 

l<M*^)l< 



1 



(t + r)(t-r) s 



l<M*,r)[< (f + r) | t _ r)3 . 



(32) 



(33) 



so we see that cftss becomes subdominant to (ft^A regarding the late time asymptotics. 
Next, in spherical symmetry □(foe = f<9_<9 + (r03c) and hence 



4>3c(t,r) 



1 

4r 



t+r 



t-r 



du / cfc d v (rf3 2 d v <pi). 

i-r J-(t+r) 



(34) 



The inner integral is an integral of a total derivative of a compactly supported function of v. The 
integration range is bigger than its support [-R, R] so the integral vanishes. This gives (ft^c = 
for t > r + R. 

The last component (ft^n = for t > r + R because the source vanishes identically in the 
integration region when inverting the wave operator as we did above. 
Finally, for late times, in the region t > r + R, 

1 



\Mt,r)\ = \Mt,r)\<- +r){t _ r)2 



(35) 



and other components of (ft can be a priori estimated to be subdominant. Therefore, in the asymp- 
totic analysis, concerned solely with the leading order behavior, it is sufficient to keep only the first 
source term in (jSJ). 

B. Higher dimensions 



Here, we want to show that the second source term in (j!4j) is asymptotically subdominant with 
respect to the first one. In dimension d+1, the estimates for 4>\i m i an d P2 and their derivatives 
can be obtained analogously to the 3+1 case and read for < t — R < r < t + R 

\Mt,r)\, \Mt,r)\, \Mt,r)\ < |^| < < -^j, (36) 



K(*,r)|<C^p \P2(t,r)\ <^r, 



\d+fo\ < ^p, 



\d-fo\ 



1 



~ {t) d ~ l 



(37) 



while all these functions vanish for r < t — R. It allows us to control both source terms in (fl4 



1 1 



< 



+ 



(t) d - 1 (r) 1 + l 



• + 



< 



r Wd-l/ r \l+l ~ / r \3+3i 



(d-i; 



< |5+/3 2 a+0i| + |/3 2 5 + 0i| + 
1 1 

< 1 + 

~ (^d[r) 2+l (t) d ~ 1 {r) 3+l r(£) d_1 (r) 2 +' ~ (r) 5+3r 



< 



(38) 



(39) 
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Both estimates are optimal (c.f. the explicit asymptotic expressions in Section IV of BCR), hence 
the second therm is indeed subdominant what justifies neglecting it in the leading order calculations. 

However, there is a problem in d+1 dimensions which is absent in 3+1. We lack an optimal 
decay estimate for the wave equation in d+1 dimensions, given a source with prescribed decay. 
The presently best known estimates (c.f. @, @]) lose I powers in the late-time decay relative to 
what is optimal. Therefore, for the above sources we are able to show rigorously only the decay 



while l/t 3+3 ' is optimal for late times. However, it does not change the fact that in the asymptotic 
analysis it is the first source term in (|14j) that dominates the asymptotics of ^3 at late times. 
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